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Abstract. Quantum lattice solitons in a system of two ultracold bosons near 
Feshbach resonance are investigated. It is shown that their binding energy, effective 
mass, and spatial width, can be manipulated varying the detuning from the Feshbach 
resonance. In the case of attractive atomic interactions, the molecule creation stabilizes 
the solitons. In the case of repulsive interactions, the molecule creation leads to the 
possibility of existence of bright solitons in some interval of detunings. Due to quantum 
fluctuations the distance between the atoms is a random quantity with the standard 
deviation larger than the mean value. 
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1. Introduction 

Many-body phenomena in ultracold quantum gases is a subject of extensive ongoing 
research. Interaction between atoms plays a crucial role in many situations and is 
responsible for the most striking experimental observations of solitons in Bose-Einstein 
condensates trapped by harmonic pQ and periodic [2] potentials. The typical number of 
atoms in such solitons varies from several hundreds to tens of thousands. Therefore these 
structures can be well described by the mean-field Gross-Pitaevsky equation. However, 
the mean-field theory does never provide an exact description of interacting quantum 
systems (e.g., due to unavoidable depletion) and it becomes interesting to investigate 
quantum effects in the soliton propagation [3]. If the atoms are loaded into the optical 
lattice, the interaction effects become more important and in the limit of the small 
number of atoms an exact quantum analysis of the solitons [HE] reveals strong deviations 
from the results provided by the Gross-Pitaevsky equation with lattice potential [S] or 
its discrete version [U U\ ■ 

The use of Feshbach resonances to control interaction between ultracold atoms 
in optical potentials is a widely spread technique allowing transformation of atoms 
into molecules and changing magnitude and sign of the effective scattering length 
of the atoms (see, e.g., [HI El)- Coherent solitons in condensed atomic- molecular 
mixtures without optical lattice were studied in several papers, see, e.g., [121 Ej. The 
mathematical model previously used for the coupled atomic-molecular condensates is 
equivalent to the model describing parametric interaction of photons in quadratically 
nonlinear crystals It was demonstrated in the mean- field limit that the resonant 
atomic-molecular interaction serves as a mechanism responsible for supporting bright 
solitons in the case of repulsive bosons and for preventing collapse in the case of 
attractive bosons |TTJ [12]. The quantum atomic-molecular solitons in the system 
without periodic potential were also studied ^QJ [E]. Optical parametric solitons in 
the system of coupled waveguides, playing the role of a periodic potential for photons, 
were recently observed experimentally [Hj. Atomic-molecular solitons in a deep optical 
lattice have been theoretically considered in the mean- field approximation JS], which 
is mathematically equivalent to the system studied in |14j . 

In this work, we demonstrate existence and study properties of the quantum atomic- 
molecular solitons in an optical lattice near Feshbach resonance. Under the quantum 
lattice soliton we understand the quantum state of the system of interacting particles 
with the localized eigenfunction and the discrete energy level belonging to a spectral 
interval forbidden for the spatially extended periodic states jlj. Note that the discrete 
energy levels belonging to the intervals forbidden for the linear waves are also a generic 
feature of the classical lattice solitons. Advances in manipulation of ultracold atomic 
systems with small number of particles per lattice site ^E] as well as in cooling and 
trapping of single atoms [T7j allow one to hope that quantum lattice solitons will soon 
become relevant for experimental research. 
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2. Hamiltonian 

We consider two atoms of mass m in an optical lattice created by a far-detuned standing 
laser wave. If the laser wavelength is Al = 2n/k L , then the lattice constant d = Al/2. 
It is convenient to represent the amplitude of the periodic potential in the form Hlj-rs, 
where u>n = hk^/ (2m) is the recoil frequency and s is a dimensionless parameter. In the 
case of a deep optical lattice every lattice site can be described by a harmonic potential 
with the frequency uj = 2u-&-\/s and the lowest-band atomic Wannier function is well 
approximated by a Gaussian with the characteristic length Z a = yjh/muj. The atoms 
in the lattice are subject to the magnetic field B, with B = Bq corresponding to the 
Feshbach resonance of the width AB. 

There are several processes which are to be taken into account in such a system: 
atomic interaction, molecule production and atomic and molecular hopping. Taking 
into account only the hopping between the nearest lattice sites as well as on-site atomic 
interactions and in the lowest-band approximation the Hamiltonian of the system is 
given by [HI d 

H = -t a ^2 a \ a 3 ~ *m 12 b t b j + ( 5 ~ 12 b l b i + ~Y 12 a \ a \ a i a i 

+ 9 12 ( b i a i a i + a i4 b i) , (!) 

i 

where a\ (b\) and a { (b { ) are creation and annihilation operators of a single 
atom (molecule) at a lattice site i, 8 — Afi(B — B ) is a detuning from the 
Feshbach resonance. Here, A/x is the difference in magnetic moments of the two 
atoms and a molecule. The atom- molecule conversion is determined by g = 
frw27rab g A.BA/i/m/(27r/a) 3/ ' 4 and the background on-site atomic interaction parameter 



is C/bg = j2/irhu (flbg/O with a bg being the background scattering length. In the 
Gaussian approximation, the atomic and molecular tunneling matrix elements are given 

by t a ,m = ^f I - (^) 2 (i^) 2e " (AL/4 ' a,m)2 - Since = U/V2, the molecular tunneling 



rate is much smaller than the atomic one. 



3. Solution of the on-site problem 



The on-site problem for the Hamiltonian {1} can be easily solved analitically. In the 
case when the atoms are on the same lattice site there are two eigenmodes which are 
superpositions of the two-atom and molecular states with the energies 



\ 



'5' - U h 



(2) 



and the probability to find a molecule 
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Figure 1. Eigenenergies in the case of two atoms on the same lattice site. Solid lines 
show the results given by Eq. which corespond to the lowest-band approximation. 
The results for the infinite number of bands [Eq. are shown by dashed lines. 
C/ bg = 0, 2^Fg 2 /(M 2 =0.1- 



where 5' — 8 — ^Tuo is an effective detuning. 

The two-atoms on-site problem was exactly solved in Ref. ^H] for the infinite 
number of bands neglecting the atom-atom interaction. The eigenenergies E are shown 
to be determined by the equation 

2y5ff» T{-E/m») 
hu T{-E/2Twu-l/2) ' 1 ; 

The eigenenergies given by Eq. (J2J) for {7b g = and Eq. (J3j) are plotted in Fig. [T] As 
we see, our lower-branch solution E_ in Eq. (j2J is in excellent agreement with the 
corresponding branch of Eq. (j3J) for arbitrary 5. The upper-branch solution E + fails to 
reproduce the second branch of Eq. P|) if 5 is far above the Feshbach resonance where 
the contribution of the second band becomes significant, remaning however in a very 
good agreement near the resonance and below it. This implies that the lowest-band 
approximation is valid if the effective detuning 5' is less than the gap between the two 
lowest Bloch bands, which is the quantity of the order of tkv, and/or if we are interested 
in the eigenmodes of the Hamiltonian (JJJ with the energies less than the energy of the 
second Bloch band. The latter is always the case in the present work. In addition, the 
parameters U\> g and g must be much smaller than the bands separation which is also 
fulfilled. 



4. Eigenmodes of the complete Hamiltonian and the soliton band 

We consider a one-dimensional model with L lattice sites and assume that L is oddf. 
Under periodic boundary conditions the eigenstates of the Hamiltonian (0) are 

m = c^Y, {Tin) 3 ' 1 1 i m o . . . o) + 4 k £ (f/n) 3 ' 1 [20... o) 

3=1 3=1 

+ 4 fc E(T/r fc )^ 1 |110...0) + ... 

3=1 

| In the case of even L there will be only unessenstial modifications in the equations. 
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+ <? { L-x)/2,k E ( f / T ^ 1 110 • • • 010 . . . 0) , (5) 

3=1 

where |l m ... 0) is a state with one molecule on the first lattice site and all the other 
sites being unoccupied, \ni . . . til) is a state with rij atoms on site i, i — 1, . . . , L. T 
is the translation operator which has the eigenvalues = exp (mk/k-i) with the wave 
number k = ki,2i>/L, v = 0, ±1, . . . , ±(L — l)/2 [3]. The eigenvalue problem for the 
Hamiltonian (^) can be written down in the matrix form 

where e™ = 5'— 2t m cos (irk/ki,). The vector A has a length (L+l)/2 and its nonvanishing 
element is A\ = \/2g. The nonvanishing elements of the tridiagonal (L + l)/2 x (L + l)/2 
matrix Qfc are given by 0] 

Qli = t/bg , Q21 = QJa = -* a V2(l + 7*) , (7) 
Qi+i,i = Qh+t = -*a(l + r fc ) , i = 2, . . . , (L - l)/2 , 

Q(L+l)/2,(L+l)/2 = -<a 



(L+l)/2 (L-l)/2- 



z 0k) ■ ■ ■ ) c (L-l)/2,fe 



, and 



The eigenvectors in Eq. (jHJ) consist of two parts c™, c^; = col 
satisfy the normalization condition 

(L-l)/2 

K\ 2 + E Kl a = 1- (8) 

i=0 

In the absence of the molecular mode, the eigenvalue problem (JHJ reduces to that 
one solved in Ref. jl], where it was shown that in the case of attractive interaction the 
energy spectrum consists always of (quasi) continuum band and a discrete level below the 
(quasi) continuum which corresponds to the bright soliton. Its characteristic feature is 
that |cQ fc | 2 ^> \c^ k \ 2 , i = 1, . . . , (L — 1)/2, i.e., the probability of finding two atoms on the 
same lattice site is much higher than all the other ones. This localization corresponds 
to the soliton solution of the discrete nonlinear Schrodinger equation and, therefore, the 
discrete level can be called a "soliton band" jlj. Our aim is to investigate the influence 
of the molecular mode on the soliton band. 

After the eigenvalue problem (jHJ) is solved, one can calculate the soliton binding 
energy which is defined as the difference of the energy at the bottom of the 
(quasi)continuum and the soliton level at v — which corresponds to k = k Q = 0. 
The effective mass m* can be worked out using a quadratic approximation for the 
eigenenergy at some small value of v (e.g., v = 1) E kl = + h 2 k\/ (2m*), which leads 
to 

m* = 2h 2 kl/[(E kl -E k0 )L 2 ] . (9) 

According to Eq. (JSJ the distance between the atoms w k is a random variable which takes 
the values Wki = 0, 1, . . . , (L — 1)/2, with the probabilities |c^| 2 / fl — |c™| 2 ). Thus, it is 
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Figure 2. Energy eigenvalues of the Hamiltonian 0J. The parameters are s = 5, 
2 v /7r<7 2 / (fiu) 2 = 0.1, S'/Jiui = 3, a,h g /\L = —0.005. Dots are the results of numerical 
solution of Eq. |JBJ for L = 41 and the solid lines correspond to the limit L — > oo. 
The spectrum is truncated from above in order to be consistent with the lowest-band 
approximation. 
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Figure 3. Probabilities of molecular and atomic states corresponding to the soliton 
band: \c'q\ (m), |cf | (i), i — 0,1,2 [|cf | is shown by the dashed line]. The 
parameters are the same as in Fig. 
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necessary to calculate not only the mean interatomic distance (w^) but also its standard 
deviation 

(L-l)/2 .j | a i2 

= - (w k )* , (w l k ) = E ( 10 ) 

i=0 J- I c fe I 

and the soliton width can be defined as yj (w%) . 

We have solved the eigenvalue problem (JHJ) numerically for finite values of L and 
analytically in the limit of infinite lattice. The results are presented below. We 
consider the cases of attractive and repulsive atomic interactions and concentrate on 
the properties of the lower-energy modes. 

4-1. Attractive atomic interaction 

We consider first bosons with attractive interactions (£/b g < 0). If 5' is negative and its 
absolute value is very large, the coupling between the molecular mode and the atomic 
mode is negligible and we have two discrete levels below the (quasi) continuum. The lower 
one corresponds to the pure molecule and another one to the atomic bright soliton. If 5' 
increases, i.e, we come closer to the Feshbach resonance, both discrete levels approach 
the (quasi) continuum. At some critical value of 8' = S'_ the upper level merges with the 
(quasi) continuum. In the numerical calculations it is not quite clear how to determine 
5'_ exactly because there are several possibilities to define it. Analytical analysis in 
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Figure 4. Soliton binding energy. The parameters are the same as in Fig. and 
k = 0. 
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Figure 5. The ratio of the soliton effective mass m* to the effective mass m at the 
bottom of (quasi)continuum band. The parameters are the same as in Fig. [2] 




Figure 6. Mean interatomic distance (wk) (left panel) and its standard deviation a w k 
(right panel) . The parameters are the same as in Fig. [3 



the case of the infinite lattice shows that the mergence occurs if at least one of the 
inequalities 

\c$ k \ > |4| ,t = 2,3,..., (11) 

is violated. We adopt this as a definition of S'_ and by doing numerical diagonalization for 
k = and for the values of parameters in the caption of Fig.Elwe obtain 5'_ = — 1.531 Uu. 
In order to have inequalities (fTT]) again fulfilled, one has to increase 8' up to S' + . Using 
the same values of the parameters we get 5' + = — 1.490 hu>. If 5' > 8' + , a discrete 
level appears above the (quasi) continuum, while the lower one which becomes a linear 
combination of atomic and molecular states remains below (see Fig- EJ) - If we increase 
5' further and go far away from the Feshbach resonance (5' 3> tkv), the contribution of 
the molecular mode into the lowest-energy eigenstate becomes negligible (FiglBl) and we 
have a pure atomic bright soliton below the (quasi) continuum jl]. The upper discrete 
level is located very far above the (quasi) continuum and cannot be interpreted within 
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Figure 7. Probabilities of molecular and atomic states corresponding to the soliton 
band: |c™| 2 (m), \cf Q \ 2 (i), i = 0,1,2 [ | c^q | 2 is shown by the dashed line]. abg/Ai = 
0.005 and the other parameters are the same as in Fig. [21 
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Figure 8. The ratio of the soliton effective mass m* to the effective mass m at the 
bottom of (quasi)continuum band. The parameters are the same as in Fig. |7| 

the lowest-band approximation. 

The soliton binding energy E\> is shown in Fig. HJ Due to the large contribution of 
the molecular mode near the resonance the binding energy is larger than its asymptotic 
value at 8' — > oo. The effective mass m* is also larger at smaller values of 8' (see Fig. EI) 
because due to the fact that t m <C i a the effective mass of the molecule is much larger 
than the atomic effective mass. The corresponding contributions of the molecular and 
atomic states into the soliton band are shown in Fig. El 

The mean interatomic distance (wk) as well as its standard deviation o~ w \. are 
shown in Fig. El The interatomic distance is well below the lattice constant d and the 
maximal localization is achieved at the edges of the Brillouin zone. However, quantum 
fluctuations are very strong and a wk > (wk). 

4-2. Repulsive atomic interaction 

In the case of repulsive atomic interaction (Ub g > 0), the situation is quite different. 
There can be only one discrete level below the (quasi) continuum which is occupied by 
the molecule as long as 8' is negative and its absolute value remains very large. Above 
the (quasi) continuum, there is another discrete level corresponding to the atomic bright 
soliton. If we increase 8' the lower discrete level approaches the (quasi) continuum and 
the probabilities of the atomic states become larger meaning that the system enters 
the bright soliton regime supported by the molecule creation (see Fig. Q). Inequalities 
(jllj) are satisfied in this regime. If we increase 8' further and reach the value 8'_, the 
probability of the molecular state becomes very small. Inequalities (jllj) are violated and 
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Figure 9. Mean interatomic distance (wf.) (left panel) and its standard deviation (j w % 
(right panel) . The parameters are the same as in Fig. 
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Figure 10. Deviations of the eigenvalues of Eq.© for finite L from that obtained in 
the limit L — > oo. A t is shown by the dashed line. The parameters are the same as in 
Fig.0 

the discrete level merges with the (quasi) continuum, i.e., the bright soliton is destroyed. 
The probability |c™| can be also interpreted as a relative population of the molecular 
component. It is a decreasing function of 5' like in the case of classical atomic- molecular 
solitons [TH] . 

If the detuning is further increased up to 8' + , the soliton band appears above the 
(quasi) continuum. For the values of parameters used in our numerical estimations, 
5'_ = 1.479 hu and 5' + = 1.521 Hcu. 

The soliton binding energy is again a decreasing function of 5' which vanishes 
at 5' = 8'_. The effective mass m* equals to the effective mass of the molecule for 
large negative 5' and reaches the value m at 5' — 5'_ (Fig. [SJ). If we come closer to 
5'_ the solitons become less localized especially at k = and the interatomic-distance 
fluctuations increase (Fig. EJ). According to our definition of the soliton width, its 
behavior is similar to that of (wk) and a w ^ shown in Fig. El 
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4-3. The limit L — > oo 

In the limit L — > oo, the wave number becomes a continuous variable. The energies of 
the continuum band are enclosed in the interval \E k \ < q k , where q k = 4t a cos (f ^) , and 
the coefficients c^ k , c™, in Eq. (JBJ) become negligibly small. Outside of the continuum 
band, the solutions have the form Cj k = a k b{ exp (if-^j) , j = 1, 2, . . . , oo. Substituting 
this ansatz into Eq. © we obtain the equation for the eigenenergy £ k = lim^oo E k 

£ 2 k = (U hg + U k f + q\ , U k = 2fj (£ k - eH . (12) 

Note that the quantity U\> g + U k plays a role of the effective atomic interaction. 
The values of a k and b k corresponding to a certain £ k are given by a k = V2(% k , 
bk = (Uhg + Uk — £ k ) /qki and the expressions for the probabilities of the state with 
two atoms on the same lattice site and the molecular state take the form 

|c \| 2 = (l-b 2 k )/[l + bl + (l-bl)S k , 

K\ 2 = 2? / (£ k - el«) 2 , (13) 



where S k = 2g 2 /(£ k - e^f . 

Eq. (|T2j) can be multiplied by (£ k — e™) 2 and treated as quartic equation for £ k 
which contains always four roots. However, depending on the values of the parameters 
only one or two roots are real and provide normalized eigenstates implying that the 
others are unphysical and should be rejected. The normalization condition (jHJ) requires 
a\/2 < 1 as well as \b k \ < 1. One can easily show that in the special case t a = t m = 
the physical solutions of Eq. (fT2j) are given by (J2J). 

We substitute £ k ± = ±q k corresponding to the edges of the continuum band into 
Eq. flHD and get 

U hg + U k ± = 0, (14) 

which leads to the identity \b k \ = 1 and as a consequence to the violation of 
inequalities (|TT|) . Eq. (|T4*|) allows to obtain the boundaries 5'_ and 5' + of the interval of 
5' within which there is only one physical solution: 

S' ± = ±q k + 2t m cos (7ik/k L ) + 2g 2 /U hg . (15) 

For the values of parameters used in the numerical diagonalization, we find S'_ = 
—1.531 hoo and 5' + = —lA79huj in the case of attractive interaction, and 5'_ = 1.479 hut 
and 5' + = 1.532 hoo in the case of repulsive interaction. The values of 5'_ are in 
perfect agreement with the results of numerical calculations for L = 41, while S' + have 
small deviations from the corresponding numerical estimations. In the special case 
ts_ = t m = 0, 5 f _ = 5' + = 5^ and Eq. (fT4^) leads to the condition 

U hg - 2~g 2 /5i = . (16) 

This is equivalent to the requirement that the effective scattering length ab g (l — 
ABAfi/5'), which appears in the mean-field theory as a result of the adiabatic 
elimination of the molecular field jH] , vanishes. The calculations presented above show 
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that in the interval of the detunings 5'_ < 5' < 5' + the effective atomic interaction is 
gradually switched from the attractive to the repulsive one. 



The probabilities 
and have the form 



'jk 



j = 1,2,..., of the atomic states in Eq.(jSJ) decrease with j 



The soliton effective mass 



hi 



(2(i-l) 



OA; I 



:i + s k 



d 2 S h 



m 



dk 2 



£k + (^bg + Uk) Sk 



fc=0/ 



2t m {Uhg + Uk) Sk 



4*2 



fc=0 



is smaller than that at the bottom of the continuum 



m 



d 2 q k 



dk 2 



k=0/ 



h 2 k 2 



(19) 



The first two moments of the interatomic-distance distribution can be shown to be 

2 



1 



bl 



1 



c. 



"OA; I 



bl 



bl)/ 



bl 



In order to demonstrate the convergence to the limit L 



Fig.HUlthe quantities A s = sup fc 



7T.(S) C 

&k ~ b k 



/huj as well as A t = sup fc 



(20) 



oo, we have plotted in 
Ejp — q k /Tluj and 



&k + Ik 



j hit) for different L, where E% 



, and are the eigenenergies 



A b = sup fc 

of Eq.© corresponding to the soliton band, the top and the bottom of the quasi- 
continuum band, respectively. A s decreases exponentially with the increase of L and 
it is very small even for low values of L. The convergence for the boundaries of the 
continuum band is slower, but the limit L — ► 00 describes quite well the results of the 
numerical diagonalization already for a few tens of the lattice sites. In addition, we have 
compared the results of the calculations obtained on the basis of numerical solution of 
the eigenvalue problem ((Hj) for L = 41 which are presented in Figs. |2Klwith that worked 
out in the limit L — ► 00 and did not find any noticeable discrepancies. This is consistent 
with the exponential decrease of A S (L). 

In the absence of the magnetic field, the molecule creation is impossible and one has 
to put g = 0, 5' = 0, t m = 0, in all the equations. In this special case, the normalizable 
solution of Eq. (JT2J) is given by = sign(Ub g )^U 2 g + , which leads to the following 

expression for the effective mass m*^ = —h 2 sign(Uh s )^JU£ g + 16t^/ (4t^). These are 
exactly the results presented in Ref. jlj. The soliton band exists again for repulsive as 
well as attractive atomic interaction, but in the case of repulsive interaction it appears 
to be a highly excited mode with the energy above the continuum band. 



5. Conclusion 



Summarizing, we have investigated quantum lattice solitons in a system of two ultracold 
bosons near the Feshbach resonance. Binding energy, effective mass, and spatial width 
of the solitons, can be manipulated varying the detuning from the Feshbach resonance. 
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In the case of attractive atomic interactions, the molecule creation stabilizes the solitons 
increasing their effective mass as well as the binding energy and decreasing the width. 
In the case of repulsive interactions, the molecule creation leads to the possibility of 
existence of bright solitons in some interval of detunings analogous to the corresponding 
classical system. The presence of quantum fluctuations leads to the fact that the 
interatomic distance is a random quantity. Its standard deviation is even larger than 
the mean value. 

The classical limit of the problem studied in the present work was considered in [T^] . 
Our results for the relative populations of the atomic and molecular components are in 
agreement with the corresponding classical results. In order to understand the transition 
from quantum to classical solitions it is necessary to perform analogous calculations for 
higher number of atoms. This can be done employing the same method as in the 
present study. However, one has to keep in mind that the dimension of the Hilbert 
space increases rapidly with the increase of the particle number and the number of 
lattice sites. 
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